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Abstract 

We study the effects of the interaction terms between the inflaton fields on the infla- 
tionary dynamics in multi-field models. With power law type potential and interactions, 
the total number of e-folds may get considerably reduced and can lead to unacceptably 
short period of inflation. Also we point out that this can place a bound on the charac- 
teristic scale of the underlying theory such as string theory. Using a simple multi-field 
chaotic inflation model from string theory, the string scale is constrained to be larger than 
the scale of grand unified theory. Allowing post-inflationary generation of perturbation 
can greatly alleviate this constrain. 
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1 Introduction 



Although inflation [H [2j [3] is considered to be able to solve many cosmological problems while 
providing the desirable initial conditions for the hot big bang universe, the practical imple- 
mentation of an inflationary scenario consistently in the context of high energy theory and 
phenomenology is still unclear |E|. One of the most profound difficulties is the realisation of 
the simple chaotic inflation with power law potential [5], which fits best with most recent cos- 
mological observations [6j [7] . This simplest possibility typically requires an initial value of the 
inflaton field far beyond the Planck scale mpi ~ 2.4 x 10 18 GeV. This super-Planckian initial 
value exposes any explicit model of chaotic inflation to uncontrollable radiative corrections and 
the inflationary predictions become not reliable. 

Fortunately, this problem can be evaded in multi-field inflation models, since the Hubble 
parameter H receives contributions from all the fields participating in inflation [HI IS]. Thus even 
with sub-Planckian field values which for the single field case will not give rise to enough number 
of e-folds J\f, the duration of inflation can be long enough to achieve all the major successes we 
expect well under the theoretical control. Moreover, there exist plenty of scalar fields in theories 
beyond the standard model of particle physics such as supersymmetry and supergravity [10J. 
Hence any inflationary model based on such a theory would naturally incorporate multiple 
number of fields 

At this point, the existence of the coupled terms between the inflaton fields is of great 
importance for the inflationary dynamics. In the single field case, the interactions with most 
non-inflaton fields are usually required to be very small to maintain the sufficiently flat inflaton 
potential. In multi-field models, however, virtually all the light scalar fields may contribute and 
there is no reason to suppress the coupled terms a priori. Also, any phenomenological constraint 
for the interaction from such a consideration should place an useful, important bound on the 
parameters of the underlying theory. This is what we would like to study in the present note: we 
will address this issue using a rudimentary model of multi-field inflation, introducing a simple 
cross term which couples the inflaton fields. The potential may arise in the context of string 
theory by breaking the shift symmetry, and the interaction terms appear with their magnitude 
being suppressed by the ultra violet cutoff scale, i.e. the string scale M s . A simple analysis on 
the inflationary dynamics can place a bound on this scale. 



2 Inflation with coupled terms 

For clarity, we first concentrate on the simplest case of the two-field chaotic inflation [T2l [T3] 
with the potential 

heading = 2 m J^ + 2 m \x 2 , (1) 

and now we include an interaction terrr|§ 

Vint = ^W, (2) 

1 The inflaton candidates are also expected to carry the standard model charges to populate the observed 
particle species, especially in the context of minimal supersymmetric standard model 

2 Note that a negative sign would induce a tachyonic instability and may lead to hybrid inflation [14] and 
even dark energy 115]. 
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where g is the coupling between and x- m addition, we take larger than m x and assume 
that rrifp and m x are not too different: if the difference is too large, the inflationary phase is 
completely dominated by and virtually x is decoupled from the dynamics, although after 
inflation it may serve as a curvaton candidate [16j [17], adding additional constraints on e.g. 
the energy scale of inflation [TBI tlH] • 

In Fig. [TJ, we present several numerical results. As can be seen, when the interaction term 
Eq. (j2J) dominates, the total number of e-folds is reduced almost by half. More importantly, 
when the coupling is strong and x tend to follow the same evolution, i.e. in the field space 
the trajectory is directly towards the minimum at the origin. Also note that for g > 5 x 10 -6 , 
J\f hardly changes. 
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Figure 1: The evolution of (thick dotted line) and x (thin solid line) versus A/". From the 
upper left plot, we have set g = 0, 8 x 1CT 7 , 2 x 10~ 6 and 5 x 10" 4 . We set the initial values 
00 — Xo — Hwpi, which without Eq. (J2J) would give Af ~ 60 from Eq. (jSJ). The bare masses 
are set to be = 10~ 5 m P1 and m x = m^/10, which gives V 1 ^ 2 ~ 5.58696 x 1CT 5 at 60 e- 
folds before the end of inflation. Note that with large number of the fields participating in the 
inflationary dynamics, significant curvature perturbation may be generated after inflation with 
the observed spectral index [20] as well as non- negligible isocurvature one [21] . The masses then 
can be quite light, but M still remains more or less the same. The total number of e-folds until 
the end of inflation is, again from the upper left plot, 62.0865, 50.3970, 39.6637 and 30.5004, 
respectively. 

We can understand the evolution of and x as follows. Given the potential Eq. flU with 
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the interaction term Eq. (j2j), the equations of motion of (f) and x are written as 



+ 3tf0+K + sY)0=O, (3) 
X + 3H X + {m 2 x + g 2 <j) 2 )x=0, (4) 



respectively, where the Hubble parameter H during inflation is given by 

1 fl 1 



H 2 = ^( ~J 2 + - n X 2 + V) , (5) 



with V being the sum of Eqs. ([I]) and (J2J). From Eqs. (jBJ and (HI), we can see that the "effective" 
masses are composed of the bare masses and the interaction, i.e. 

^J(eff) =™ 2 + #V , (6) 
m l(es) = m l + 9 2 4> 2 , (7) 



for and Xi respectively. Thus two different effects are competing, and the dynamics is deter- 
mined by the term which is important: if the bare masses dominate, <fi and x are practically 
decoupled and it is well known that the total number of e- folds Af, neglecting 0(1) corrections, 
is given by [121 13] 

Af = [ Hdt=^±>£. (8) 

Conversely, when the interaction is strong enough, this prediction becomes different: the equa- 
tions of motion, especially at the early stage of inflation, are simplified to, with the replacement 

P 2 = 4> 2 + X 2 , 

p + 3Hp + g 2 p 3 = 0, (9) 

i.e. the equation becomes that of the quartic potential and both and x essentially follow the 
same evolution. In this case, M is given by 

^ = 4^, (10) 

8mp[ 

and is reduced by half compared with the case where the bare masses are dominating, which is 
clear from Fig. [lj This is also why the coupling larger than a specific value does not alter M 
appreciably: this value corresponds to what makes the interaction term larger than the bare 
mass contribution. 

The case of intermediate interacting energy should be of most interest. Since the slow roll 
approximation is valid for most duration, we can write the equations of motion for <fi and Xi m 
terms of Af, as 

ti + rr&Q=0, (11) 

where we have used H 2 w V/(3mpi), a prime denotes a derivative with respect to A/" and the 
subscript i stands for (p and %■ We can thus write the solution of Eq. (ITTj) as 



-v^AAmpi, (12) 
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where the slow-roll parameter 



mpi (V, 



2 \V (13) 

is assumed to be nearly constant, which is a good enough approximation in the slow roll regime 
of inflation. Since we are assuming that the energy associated with the heavier field, <fi here, 
and the interaction energy are relatively large, the number of e-folds until 4> contributes follows 
that of a quadratic potential of 4> which is independent of the effective mass of (ft given by 
Eq. (|SJ), and is estimated simply as 

**~A- (i4) 

Until then the lighter field, Xi evolves according to Eq. (fl2|) . The amplitude of x a t the moment 
when exits inflationary regime, x<j>i is thus given by 



^-Xo-^-r^— • 15 
4 mpi 

Afterwards only x drives inflation. Hence the number of e-folds J\f x , achieved during this phase, 
is simply given by the standard result 

K « -^r > ( 16 ) 

where we have assumed for simplicity that inflation ends when % = 0. The total number of 
e-folds is then, from Eqs. ( jMl) and (|T6|) . 



Af « N<j, + M x 
1 



4m P1 



(17) 



For example, let us take the numbers from the upper right plot of Fig. [TJ With the given 
numbers, roughly 



1 

1 



-ml4>l ~ 10 9 m P1 , 



2„,2 m-H™4 



PI ' 



i^xg ~ io- 9 ™ 4 P1 • (18) 

Thus, from Eq. (fill . A/^ ~ 30. Meanwhile, using Eq. ({TBI , we can estimate 

= (m 2 x + 9 2 <fi)Xo 

x vx ™\$l + m \xl + a 2 (t>lxl 

_ 2 Xo™pi 
~ 9 ~ 

m 4> 

^0.07, (19) 
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where we have used our assumption that the potential energy of <p and the interaction energy 
are of the same order of magnitude and are very large compared with the potential energy of 
X, as can be read from Eq. (jTSl) . Substituting Eq. (Tl9l into Eq. (ITTI) . then we have 



J\f & 50. 



(20) 



which is in good agreement with the numerical estimate Af ~ 50.3970. Also we note that both 
the bare masses and the coupling do not affect the duration of the inflationary epoch but they 
determine the amplitude of the power spectrum V, see Eqs. (l26l) and (|28l) . In Fig. we show 
the dependence of Af and V on the coupling g. 
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Figure 2: The plots of (left) Af and (right) V 1 ^ 2 versus g. Square dots denote numerical results, 
and solid lines are interpolated curves. We have used the same masses given in Fig. (TJ but 
this time the initial values are set to be <fo = Xo = 16mpi so that even when the interaction 
completely dominates the potential, with Af given by Eq. (flUj) . we still have more than 60 
e-folds. In the right panel, we have evaluated V 1 ^ 2 at 60 e-folds before the end of inflation. 
As can be seen in the left panel, Af is reduced almost by half when the interaction becomes 
strong. Note that for a wide range of g, Af remains more or less the same and only the value 
of g taken from a limited range can significantly modify Af. The constraint on g beyond this 
range can be obtained from the magnitude of V 1 ^ 2 , which is approximately 5 x 10~ 5 by recent 
observations [3[7]. V 1 / 2 is linearly dependent on g in logarithmic scale, which suggests that 
Eq. ( l28l) is indeed the case. According to the above numerical results, g ~ 10 -6 gives the 
observed value of P 1 / 2 with the effect of the interaction being prominent but not completely 
dominating: see the lower left plot of Fig. [U We have also calculated the spectral index n s , but 
irrespective of the value of g it remains between 0.95 and 0.97. From Eqs. (127]) and (1291 . it is 
expected that n s is almost independent of g. 

We can easily expand the argument to include arbitrary number of fields: there are total N 
fields (pi, each with the mass so that the leading piece is 



N 



leading 



(21) 
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with the interaction between fa and fa (i ^ j) being given by 



Va = \g\A$ • (22) 



The potential is then written as 



N 1 1 

^ = E^ 2 + Ek^?- ( 2 3) 



2 2' 

i=i 

Then the equation of motion of fa is given by 

fa + 3Hfa + I m 2 + 9%$ ) 0, = . (24) 

Thus while the bare mass remains fixed, the number of interactions can become very large 
according to the total number of fields, making them completely dominating unless the couplings 
are all exceptionally small. For simplicity, let us assume that all the masses and the coupling 
constants are the same, and write them as m and g, respectively. Then with N ^> 1 so that iV 
and N — 1 are more or less the same, Eq. ff24|) is reduced to 



p + 3Hp + [m 2 + Ng*p 2 ) p = . (25) 

Compared with Eq. (Q, when the interactions are strong, g should be suppressed by a factor 
of v^/V to maintain the same predictions. 

Now we briefly mention perturbations. We may follow the SAf formalism [22l [23l EH [25] , or 
the standard procedure for the single field case j26j E7] , to calculate the observable quantities 
such as the power spectrum V and the spectral index n s : if the masses are dominating, we have 

2 4 

P= T-iTi ( 26 ) 



= (27) 
P 2 



or when the interactions are important, 



,2, 



V = 7^2 6 » (28) 
76o7T TTlpj 

24m?,, , , 

n s - 1 = ^ . 29 

When the masses or the couplings are different, typically we have a redder spectrum [28J. At 
this point, it is interesting to note that the results are confined between those of the fa 2 and </> 4 
theories. This can be seen from the geometry of the field space: the interaction gf^fajfa 2 lifts the 
whole potential, with the steepest rise along the line fa = fa which corresponds to the quartic 
potential, while leaves fa = and fa = 0, i.e. the mass eigenstates, intact. Thus in the field 
space the descent is mildest along each axis with quadratic dependence, and is steepest along 



6 



4>i = 02 = ■ ■ ■ = 4>n with quartic dependence. Thus at the early times the trajectory tends 
to be towards the origin along the steepest descent. This geometric consideration should be 
applicable to generic power law potential and interactions. Also we should note that the shape 
of V is saw-edged. With not equal masses (or couplings in general) heavier fields will drop 
out of the inflationary regime first. Near the moment of a single drop out, we can expand the 
potential then a change of slope is experienced: with A denoting the slope, when iV-th field is 
relaxed to minimum and quits the inflationary regime, we have J29J [30] 



AA 



A 



m 2 N (j) N 



(30) 



where the right hand side is evaluated at the moment of the drop out. This will lead to an 
enhancement of V across the corresponding scale, which becomes mild when the masses are 
densely spaced. In Fig. [3] we plot the phase portrait of x of Fig. [TJ from which this momentary 
feature can be read. It is also noticeable that with strong enough interactions, the level of 
non-Gaussianity in general will not be significant, /nl during inflation is known to become 
large when the trajectory in the field space is highly non-trivial [31]. For the case where the 
interaction terms are dominating, however, all the fields evolve essentially in the same, or at 
least similar way, making the trajectory in the field space varying very smoothly. This will 
lead to mild, even negligible enhancement of /nl- Although it should be possible to obtain 
considerable non-Gaussianity with a very special choice of the initial conditions, but this looks 
not quite likely. 
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Figure 3: The evolution of x m the phase space. The values of the parameters are taken from 
Fig. [T]so that the cases with g — (black), 8 x 10~ 7 (green), 2 x 10 -6 (blue) and 5 x 10~ 4 (red) 
are shown. It can be seen that x tends to follow an attracting trajectory and that thus the 
slow roll approximation is good enough during most duration of evolution. The 'kinks' of the 
green and blue lines denote the moment of drop out of <ft, so that from then x quickly follows 
the uncoupled evolution. 

Finally, we mention the possibility of parametric resonance due to the oscillation of a field 
near its minimum. In the theory of preheating [321 13~3] . the strength of parametric resonance is 
determined by the parameter 

o 2 $ 2 

9 = kr (31) 
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where $ is the amplitude of oscillation of the inflaton <\> when the inflationary phase ends. 
Only when q ^> 1, the field coupled to 0, say X-, g° es through many instability bands of 
the Mathieu equation and the resonance becomes prominent. In the present case, however, 
because of the relatively large Hubble parameter, the energy of the massive fields which would 
oscillate earlier in the absence of other light fields is dissipated and those fields are overdamped, 
showing negligible oscillation. Even for the light fields which oscillate later, it is known that the 
amplitude near the end of inflation is suppressed by a factor proportional to y/~N [31], making 
the amplitudes of oscillation far smaller than that of the single field case. Hence we have small 
q and the resonance is suppressed. If we increase g to make q large, as discussed before all 
the fields follow the same evolution and oscillate in phase, heavily suppressing the possibility 
of parametric resonanc^f]. We have tried different parameter ranges but have not found any 
instability of x- This seems to suggest that we need fine tunings to observe any parametric 
resonance. 



3 Bound on string scale 

Now we turn to our next point how the discussion in the previous section may put a bound on 
the string scale M s . It was argued [36] that in string theory the ubiquitous string axion fields 
can give rise to the uncoupled quadratic potential, Eq. (pi) , as the leading approximation, 
provided that the fields are displaced not too far from their minima. The potential is a sum of 
the axion contributions^, 



heading = A i 1 ~ COS \ 



N 

l •>—, i \ 

(32) 



where A 4 is the axion potential scale and fi is the axion decay constant. Multi-instanton 
corrections generate the coupling terms between the axions, 

where M is the cutoff scale of the theory, which in string theory should correspond to the string 
scale M s . Typically Aj is smaller than mpi by many orders of magnitude thus provided that 
Aj <C M s Eq. (|33|) is highly suppressed and can be ignored so that Eq. (|2~T!) is a good enough 
approximation of Eq. (|32|) . 

For small field values, we can expand Eq. (1321) and easily find that 

m = , (34) 

Ji 



3 There is, however, some possibility in the context of A</> 4 theory, which is equivalent to the present sit- 
uation [35j . However pursuing this idea is beyond the scope of the present paper, and we do not discuss it 
here. 

4 The next higher order self-interaction near the minimum will then be negative, — (2wAi/ fi) 4 c/>f /24. Thus we 
do not consider these contributions. But even if we take them into account, still the argument here is applicable. 
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and from Eqs. (1221) and (1531). 



m,m,- r- 9 A- A- 

'"^vkrr^im- (35) 

Moreover, it is known [37] that in string theory is likely to lie within 

0</i<m P1 , (36) 

so that any field value beyond m-p\ seems prohibited. But nevertheless as mentioned in the 
previous section with a large number of fields it is possible to have the simplest chaotic inflation 
with quadratic potential as Eq. ( 12~Ti) . which is favoured by recent cosmological observations. The 
important point is that, still the predictions generally follow those of the single field model, so 
that for example the overall mass scale is constrained to be 

(rn) ~ O (l(T 5 ) topi (37) 

from the observed amplitude of density perturbations [38] • This in turn can put a reasonable 
bound on the string scale M s . 

For simplicity, let us assume that all the masses and the couplings are the same and that 
all the axion decay constants are set to be mpi. The observed density perturbations require 
Eq. (1371) . thus using Eq. (IMjl we have 

m = Kr 5 m P1 , 38 

mpi 

so that the axion scale is 

A ~ Kr 3 m P1 . (39) 

Now, using Eq. fl35l) . we can estimate 

>-.L() /'^P1 X 



[lil) • (40) 

thus an appropriate bound on the couplings can constrain the possible range of the string scale 
M s : for example, from Fig. [T]let us take the bound g < 10~ 6 for the couplings not to completely 
dominate the dynamics. Then , we obtain 

M s > 10~ 2 m P i ~ 10 16 GeV ~ M GU t • (41) 

Thus, in this case the string scale smaller than the scale of grand unified theory Mgut would 
lead to comparably short period of inflation so that it may not solve the cosmological problems. 
In this context, M s as large as m P1 highly suppresses the couplings and is thus favoured. 

We can obtain the same estimate by comparing the magnitudes of Eqs. (I32p and (1331) as 
follows |39j . By assuming that all the A, and fa are of the same order of magnitude for simplicity, 
we have 

heading ~iVA 4 , (42) 

V im j^- , (43) 
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where we have used the fact that for a random distribution of 9 between and 2tt, (cos 2 9) = 
1/2 ~ 0(1). The interaction energy Vj nt increases proportional to iV 2 since it is summed over 
two different indices. Since we expect that inflation is driven dominantly by Vj ea dmg, which 
should be the energy scale of inflation, we can write 

heading M g 

V int 

That is, the inflationary energy scale is bounded from above by the cutoff energy scale of 
the underlying theory, which is quite reasonable. To reproduce the predictions of the chaotic 
inflation with quadratic potential, with the inflaton mass given by Eq. (1371) . we have M- m f ~ 
Mgut about 60 e-folds before the end of inflation. This again gives M s > Mgut- In Fig. @J we 
show the dependence of Af and V on the string scale M s . 



M, 



M 



inf 



> 1. 



(44) 
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Figure 4: The plots of (left) J\f and (right) V 1 ^ 2 versus M s /m-p\ with the same potential and 
parameter values. As in Fig. [2j square dots are numerical results, and solid lines gives the 
interpolated curves. 

An important point we can read from Fig. [4] is that, the coupling g cannot be indefinitely 
small. This is because g is inversely proportional to M 2 as Eq. (135|) and the current ultra violet 
cutoff scale is bounded by m-p\ from above. We immediately face a difficulty regarding the 
amplitude of V, which becomes unacceptably large as M s becomes even slightly less than mpi 
provided that the masses are of the right magnitude which gives "P 1 / 2 ~ 5 x 10~ 5 , i.e. Eq. (I37j) . 
In this case, as can be seen in the right panel of Fig. HI even the largest conceivable string scale 
M s ~ 0(mpi) is only marginally acceptable. 

This estimate is, however, completely changed if we allow a post-inflationary generation of 
perturbations such as the curvaton mechanism [TgllTTj. Eq. f[H7|) assumes the simplest possibility 
that no curvature perturbation is generated during the post-inflationary dynamics and that all 
the observed inhomogeneities in the universe are purely due to the fluctuations of the inflaton 
fields. If we do not require this standard lore we can liberate (m) from the constraint Eq. (|37l) 
and it can be reduced by large orders of magnitude. From Eq. (|35l) this gives rise to a large 
allowed range of M s . In that sense, a post-inflationary generation of perturbation is not just an 
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option, but is rather strongly favoured in the light of richer string phenomenology which is due 
to considerably smaller M s than mp\. If additional curvature perturbation is generated via the 
curvaton mechanism, it is required [HJ H] that during inflation a > 10 _6 mpi, or equivalently, 

H"^M^> i -^ m 2 (45) 

so that we are allowed have an intermediate inflationary energy scale of Mj n f > 10 n GeV. 
Eq. fl37|) then can be lowered as much as 

(m) ~ O (10 3 ) TeV . (46) 

Then, from Eq. (J35l) . 

3 >10-(^) 2 , (47) 

which now allows M s much lower than Mqut- This is perfectly compatible with phenomeno- 
logically interesting intermediate string scales of M s ~ 10 10_12 GeV, which may naturally arise 
in specific string compactifications such as large volume compactification [40]. From Eq. (I3~4l) 
this is reduced to the problem of constructing very large fi, which is still unclear in the context 
of string theory [37] , although there exist some possibilities [HI 02] . 



4 Conclusions 

We have studied the dynamics of multi-field inflation in the presence of the interactions between 
the inflaton fields. In single field models the inflaton field is required to be coupled very weakly 
to other fields to maintain the flatness of its potential, otherwise the successful predictions of 
the slow-roll conditions which are consistent with the most recent observations are apt to get 
spoiled. In multi-field models, in general all the light scalar fields contribute to the inflationary 
epoch and there should exist coupled terms between the inflaton fields, which presumably have 
profound effects on the inflationary dynamics. As the simplest possibility we have taken the 
multi-field chaotic inflation model as the example with quadratic coupling terms between the 
inflaton fields. As shown in Fig. [TJ the first observation is that the total number of e-folds M 
can be considerably reduced. This happens when the coupled terms dominate over the bare 
masses, since the quadratic couplings act as the quartic potential in the equations of motion 
of each inflaton field. We can follow the same reasoning for any general power law potential 
and interactions, and the observable quantities such as the power spectrum V and the spectral 
index n s vary between the predictions of the theory with the lowest power and those of the 
highest power of the potential. 

This consideration can place a bound on the characteristic scale of the underlying physics 
such as the string scale m s for string theory. We have considered the case where the ubiquitous 
string axion fields can play the role of the multiple inflaton fields. With the assumption that 
they are displaced not too far from the minimum they can reproduce the uncoupled quadratic 
potential, Eq. ([21]) . as the leading approximation, but via multi-instanton corrections different 
axion fields are coupled as Eq. (1331) . These coupled terms are suppressed by the string scale, 
thus the strength of the coupled terms is dependent on the magnitude of the string scale. Hence 
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from the requirement that the interactions be not too strong to hinder a long enough period of 
inflation, we can place an interesting bound on the string scale. We have considered here the 
simplest quadratic potential, and in this case a large M s of at least (9(Mgut) is preferred. But 
allowing post-inflationary generation of perturbation can give much more relaxed constraint. 
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